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I The zero-point energy of atoms of two stable isotopes of helium - *He and ^He - is high enough to 

prevent their solidification even at extremely low temperatures without application of external pressure. 
(~| So they belong to the class of quantum fluids with strong correlations between atoms but their behavior 

Q is quite different at low temperatures. The first one represents a Bose-system and shows superfluid 

transition at 2.17 K while atoms of 3He are governed by Fermi-statistics and superfluid transition can 
' ' be observed only at much lower temperatures (about 1 mK) when the pairing of two atoms occurs. 

But even at helium temperatures (1.5-4.2 K) the effects of quantum statistics for ^He atoms becomes 
pronounceable especially in nanoscale confinement (nanoporous media, thin adsorbed layers on solid 
substrates) and in the presence of nanoscale disorder induced for example by silica aerogel strands. In 
recent years the problem of correct description of quantum fluids in the confined geometry at nanoscale 
length has emerged. It has been recognized that the quantum fluids at these circumstances can be 
considered as a new state of quantum matter due to close values between characteristic lengths for 
these quantum liquids and the size of geometrical confinement and significant contribution from the 
surface atoms. So one has to apply new physics to describe such systems with taking into account 
their complex nature. For example, last two years the attempts to develop the fractionalized two-fluid 
hydrodynamics for nanoporous media with fractal dimensions have been made. The actuality of such 
new hydrodynamics becomes very clear for the last development in chemical synthesis of different 
• • kind of aerogels with nanopore structure as well as numerous studies of nanoporous substances. One 

^ of the interesting obtained results that density waves (the first sound) and temperature waves (the 

second sound) become strong coupled even in the absence of viscosity, so it is purely geometric effect 
of fractal space of nanopores. In the present report we will review the procedure, results and discuss 
^ the issues for this approach. 
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1. Introduction 

The zero-point energy of atoms of two stable isotopes of helium - "^He and ^He - is high 
enough to prevent their solidification even at extremely low temperatures without applica- 
tion of external pressure. So they belong to the class of quantum fluids with strong corre- 
lations between atoms but their behavior is quite different at low temperatures. The first 
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one represents a Bose-system and shows superfluid transition at 2.17 K while atoms of 3He 
are governed by Fermi-statistics and superfluid transition can be observed only at much 
lower temperatures (about 1 mK) when the pairing of two atoms occurs. But even at he- 
lium temperatures (1.5-4.2 K) the effects of quantum statistics for ^He atoms becomes pro- 
nounceable especially in nanoscale confinement (nanoporous media, thin adsorbed layers 
on solid substrates) and in the presence of nanoscale disorder induced for example by silica 
aerogel strands. In recent years the problem of correct description of quantum fluids in the 
confined geometry at nanoscale length has emerged It has been recognized that the 
quantum fluids at these circumstances can be considered as a new state of quantum matter 
due to close values between characteristic lengths for these quantum Uquids and the size of 
geometrical confinement and significant contribution from the surface atoms. So one has to 
apply new physics to describe such systems with taking into account their complex nature. 
For example, last two years the attempts to develop the fractionalized two-fluid hydrody- 
namics for nanoporous media with fractal dimensions have been mad£ill The actuality of 
such new hydrodynamics becomes very clear for the last development in chemical syn- 
thesis of different kind of aerogels with nanopore structure as well as numerous studies 
of nanoporous substances. One of the interesting obtained results that density waves (the 
first sound) and temperature waves (the second sound) become strong coupled even in the 
absence of viscosity, so it is purely geometric effect of fractal space of nanopores. 

In the present report we will review the procedure, results and discuss the issues for this 
approach. The rest of this paper is organized as follow: In Section|2]we make a little review 
of existed model of superfluid bulk helium-4. In the next Section [3] propeties and pecu- 
liarities of aerogel are discussed. In Section|4]possible fractionalization of hydrodynamical 
model is shown. In Section|5]we speculate about another way of considering complex frac- 
tal structure of aerogel for description of superfluid helium behavior In Section|6]we make 
a summarization and some novel results are given. 



2. Existed model 

There are several approaches to describe the behavior of bulk superfluid helium-4. For 
example, two-fluid modeGE, the microscopic description based on the Gross-Pitaevskii 
equatioiEHniand others. 



2.1. Two-fluid model 

The most famous model is two-fluid model (TFM) proposed by Landau in the past century. 
In this model superfluid helium with density p is considered as a two component system: 
an uncondensed, normal component with density p„ with velocity v„ and a condensed, su- 
perfluid component characterized by density — p — p„ with velocity v^ . Without dissipa- 
tive terms one finds the following system of hydrodynamical equations (so-called Landau- 
Khalatnikov equationflU ) 
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^ + div{p„\„ + p,.v, ) = 0, (1) 

d , , d 

3- (PnV„ + P.V,); + 3— n,,t = 0, (2) 

at oxk 



dS 
dt 



divS\„ =0, (3) 



d\ 

+ m(\,V)\s + VjU = 0, (4) 

ot 

The entropy is denoted by S, the mass of an atom by m, the pressure by p and the chemical 
potential by jU. In (j2]) the stress tensor Ylji^ is given by 

^ik = PnV„iV„k + PsV.siV,k + pdik, (5) 

and the sum over the index k is assumed. 

One can obtain from this equations that there are two type of collective motion: motion 
of the fluid where both components move in phase is called first (ordinary) sound, while 
second sound is associated with out of phase motion of the two components. The above 
TFM equations describe the flow properties of superfluid '^He in the bulk of the liquid. 
They are not valid in the immediate proximity of walls and free surfaces, where the effects 
of so-called "healing" are encountered. Moreover, Landau-Khalatnikov equations cannot 
be used in the vicinity of X point, where large variations of superfluid density in space as 
well as in time may occur. 

For taking into account such "healing" effects, which has quantum nature, one can 
consider some models derived from nonlinear Schrodinger equations. 



2.2. Nonlinear Schrodinger equation 

The Schrodinger equation itself is applicable only si T — OK and small interaction be- 
tween helium atoms, when all of them are in condensed state. In that case we can introduce 
macroscopic wave function of condensate. Than velocity field is determined as probabil- 
ity flow of this wave function and it is proportional to phase gradient. He-II is dense fluid 
of strongly interacted bosons and applicability of NLS is questionable, but it helps us to 
elucidate basic peculiarities of superfluid dynamics. 

If one writes short range pair interaction between particles in form U {x~x') — Uq5{x — 
xf), No = J \^\^d^x - number of particles in condensate and Vq - total condensate volume, 
than one can write nonlinear Schrodinger equatioi^in form of 

/— = — V^^+I^p-^. (6) 
ot 2 

Here average density po = A^o /^o is set to unity, unity of length is h/ -/pj'o and unity 
of time h/{poUo). 
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After substitution ^(x,f) = p{x,t)e'^^-'^-'^ and disparting real and imaginary parts in 
(|6| one can obtain 

|^+V.pv = (7) 

where v = V9. Except for high-order derivatives in right hand side of equation ([8]l, which 
can be omitted in hydrodynamical limit, system of equations (|7]i-(|8]l is equivalent to Euler 
equation for nonrotating ideal fluid with pressure defined as p{p) = p^/2. Pressure depends 
only on density p because fluid is considered at absolute zero. 

For considering dynamics of He-11 at nonzero temperature. Hills and Roberts have im- 
proved by introducing additional terms, that is proportional to the gradient of su- 
perfluid component. In that case equation (j4]) can be rewritten as 

•^v. , V7 /^v? , ^ f^^^ 2. 1 dri .^^ ,2 



where tj(Pj) - some undefined function; stress tensor has to be modified and rewritten as 
Ylij = p„v„iv„j + p, v,,v,y + T] V,psVyp, + (p- r\Ps'^^~Ps (^Ps)^^ (10) 

New terms that include spatial derivatives of p^ are responsible for macroscopic quantum 
effects, like healing lengtlPSl 



2.3. Correlated density matrix theory 

Another way for describing quantum strongly correlated system (superfluid helium is just 
one example of such systems) is correlated density matri^^ theory that provides the 
method of choice to analyze the microscopic structure of strongly correlated quantum flu- 
ids in thermal equilibrium. The main point of this theory is concept of renormalized bosons 
and fermions. Their mass can depend fully on temperature, density and wave vector So 
one may discuss the properties of strongly correlated quantum fluids at finite temperatures 
with reference to the background gas of renormalized free particles. 



2.4. Extended irreversible thermodynamics 

Recently the extended irreversible thermodynamics has been applied to describe He-I and 
He-II dnil It is shown that the behavior of helium II can be described by means of an 
extended thermodynamic theory where four fields, namely density, temperature, velocity, 
and heat flux are involved as independent fields. This model is able to explain the propaga- 
tion of the two sounds that are typical of helium II, and the attenuation calculated for such 
sounds is in agreement with the experiment results. After some approximations, this model 
is reduced to TFM, so it can be considered as some generalization of TFM. 
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3. Aerogel 

Silica aerogels are synthesized via a sol-gel process and hypercritical drying which enable 
production of tenuous solids with porosity (p as large as 99.8 % and unique acoustic proper- 
ties. Silica aerogel are known to be good examples of fractal materials. A revealed by small- 
angle x-ray-scattering (SAXS) experiments or small-angle neutron-scattering (SANS) ex- 
periments, they are made of a disorder, but homogeneous, array of connected fractal clus- 
ters resulting from the aggregation of primary particles. The analysis of the wave-vector 
dependence of the scattering intensity I{q) has permitted the determination of two charac- 
teristic length scales which are the average size aw 10 A of the particles and the average 
size of the clusters £, w 100 A. At length scales from a to ^ silica aerogels show a fractal 
behavioJSl. 

The computational confirmation for cluster structure of aerogel has been obtained by 
modeling as well as by the geometrical analysis of the diffusion limited cluster-cluster 
aggregatiorHmil 

Also it has been demonstrated^ that it is long-correlated structure of aerogel that makes 
an essential influence on liquid ^He behavior near A-point. 

3.1. Behavior of^He inside aerogel 

There has been considerable interest in the behavior of superfluid ^He in the presence 
of a random disorder induced by highly open porous media, like aerogel. Understanding 
the results of acoustic experiments is important when dealing with porous media. Use of 
liquid ^He offers unique advantages due to the existence of the superfluid phase with more 
than one sound mode. In a porous media where the normal component is clamped by its 
viscosity and only the superfluid component can move, fourth sound (relative motion of 
the superfluid and normal fluids) propagates and can be used to determine the superfluid 
fraction. 

The high-porosity aerogels are so soft that the aerogel matrix and the clamped normal 
fluid moves as the results of pressure and temperature gradients, unlike other porous media. 
This results in sound modes intermediate between first and fourth sound and a second- 
sound-like mode. In that case, proposed by Biot theory of acoustic propagation in porous, 
fluid filled, macroscopic ally homogeneous and isotropic media, is not applicable. 

It is very interesting to study the possible influence of geometrical confinement with 
fractal dimensionality on the flow properties of superfluid '^He in the framework of TFM. 
In order to shed a light on geometrical factor itself one can neglect here by quantum healing 
and any dissipative processes. 

3.2. Nonextensivity 

So hereafter we suppose that effectively aerogel can be considered as a cluster with a fractal 
mass dimension'^ the nanopores of which are filled in by liquid helium. 

Note that thermodynamic limit conditions are violated for helium atoms inside 
nanopores because of a huge inner pore surface of aerogel (up to 3000 m^/g for an aero- 
gel with density 2 mg/cm^). Namely the ratio of total number of helium atoms to the 
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total cluster volume V is not constant eA N ,V ^ °° and nonextensivity of physical prop- 
erties takes place for helium atoms in two nanopores. In this case the methods of non- 
extensive thermodynamic^ should be applied to construct the two-fluid hydrodynamic 
model and non-extensive entropy like Tsallis entropy should be introduced. We have to 
note here that beyond thermodynamic limit even the usual Boltzmann-Gibbs entropy be- 
comes non-additive, but the additivity is restored when thermodynamic limit conditions are 
taken ESl. 



4. Fractionalized and nonextensive hydrodynamical model 

In general case a macroscopic quantity Q{A,B) associated with the total system may be 
expressed in terms of the same quantity associated with the subsystems, Q{A) and Q{B)^ 

Q{A,B)=h^[Q{A),Q{B)], (11) 

where fx^ is a symmetric bivariate function depending on a constant Ag. Of course, for 
given quantity Q there exist many functions which satisfy the composability property ( fTT) . 
However, additional assumptions drastically reduce their number. For example, the ther- 



modynamic equilibrium may be used as a constraint on the form of in Eq. (111. For 
energy Ex^"^ and entropy of helium inside nanopores it leads to 

Ex{A + B)=Ex{A)+Ex{B)+XEx{A)Ex{B), (12) 
S, {A + B)^S,{A)+ S, {B) + qS, {A)S, {B), (1 3) 

where X and q are parameters of nonextensivity and are determined by the properties of 
system. Suppose that a local equilibrium between liquid ^He in different nanopores takes 
place. So we can re-define such thermodynamics quantities as temperature and pressure in 
form off! 

_ l+qS, (dEx\ 

Further, it is possible to introduce the spatial pressure and density distributions in fractal 
cluster as it has been made in our previous worlP 

Pphys{r) = Pf{r)Xp{r), (16) 
Pphys{r) ^ Pf{r)Xp{r), (17) 

where r is distance from center of fractal cluster, Pf{r) is distribution of pressure in pore, 
Xp{f) is a fractal factor-function. In the case of Euclidean space with D = 3 this factor- 
function should be equal to unit: Xpi'') = 1- Such kind of fractionalization procedure can 
be applied to any thermodynamic quantity A i.e. 

A=Af{X,Y,Z,...)xA{r), (18) 
limo^3ZA = l. (19) 
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From the defined thermodynamic quantities ( 14|15 1 and eqn. (16i the fractal factor- 
functions for energy, entropy and temperature are derived 



%E = XS 
XT = 



Xp 



l+XEf{l-Xp)' 

i+qSfXs i + XEf 

l+qSf l+XEfXE' 



(20) 
(21) 



Because all thermodynamic quantities should be expressed in terms of physical (observ- 
able) variables, one can propose the following definition of the generalized free energy 



I dE 

\dSphys 



Jphys- 



(22) 



which represents no more than the Legendre transformation. So the entropy differential 
equals to 



_ 1 +IE^ 
dSphys = - —dSq, 
l+qSg 

In the first order with respect to XEf and qSg the entropy can be written as; 



(23) 



Sphys ^ - \n{\ + qSy) + A 



1 +qSq 

Finally the fractal factor-function for Sphys is 



EliSM^^^_^^^S^^^ (24) 



Xsp — Xs ^ XpiXi 



I XEf - qSf 



XH 



(25) 



where H = 1 /Sf J EfdSf. After substitution the fractionalized thermodynamical quanti- 
ties into TFM one can derive the main system of equations for the fractionalized twofluid 
hydrodynamic model 



^Xp + diviPsfXp-Vs+Pn/Xp-Vn) = 0, 
dpfOf 

—^Xsp+divipfG/Xspyn) = 0, 

p^fXp^ = -'^^ipfXp)+psf(yfXpXa^iTfXT), 



PnfXp 



dv„ _ Pnf 



dt 



pf 



"^(PfXp) - PsfOfXpXcy^iTfXT): 



(26) 
(27) 
(28) 

(29) 



where = 1 + (z,, - 1) 



''^f +XH]. From eqns. (|26 



27 



28 



29 1 two equations for 
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Fig. 1. Profile for pressure wave (black line) and the temperature wave induced by it(dashed line). 



waves of pressure and temperature follow 

^^=Vp/ + 2 — V,,+ _,,, (30) 

\^^V^Tf{l + {x„-l )M) + VTf (mVXp -H-M)^' 



+ Tf{qSf-XEf) hh-^^ 

\ AP 
«2 Pf \^^fJpf\ Xp J 

where Mj = {dpf/dpf)sp u\ = PsfCrj / Pnf{dTf /da/jp^ are the squared first and second 
sound velocities respectively and M — {^Ef + qSf)/2 + XH, N — (3A£'/ — qSf) /2 + XH. 
On figure [T] the profiles for pressure wave and for the temperature wave induced by it 
are shown. It is seen from eqn. ( [3T| i that the coupling between pressure and temperature 
waves appears even in the absence of ^He viscosity and aerogel sceleton inertia, which is 
undoubtedly the effect of not only fractional dimensionality of nanopore space, but also of 
the nonextensive nature of thermodynamical quantities for He-II inside nanopores. 



5. Fractional nonlinear Schrodinger equation 

Let us assume that confinement (aerogel strands) forbid realization of some Feynman's 
pathes for particular helium atom in contrast to helium atom in the free space. For tak- 
ing into account such peculiarities one can propose to make use of fractional Schrodinger 
equatiorP31 

Fractional Schrodinger equation (FSE) is derived from path integration over Levy's 
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pathes and is given by 

in^^!^^Ka{-h^vY'v{r,t)+V{r,t)xif{r,t). (32) 
Here Riesz fractional derivative is introduced as 

i-h^vY'vir^t) - I dV-'/'lplXp-O, (33) 

(p{p,t) - / fl're-'>/V(r,f)- (34) 



5.1. Galilean noninvariance of fractional Schrodinger equation 



Let us consider tranformation of equation p2] i when inertial reference frame is changed. 
With such transformation spatial coordinates and time are changed as f' = f , r' = r — vt, 
where v is some constant relative velocity. Then time and fractional spatial derivative is 



rewritten aJ^Sl 



dj ={d,,-vd,.)f, (35) 

(V.2j = (V,,2j . (36) 
Let us assume that in old reference frame K FSE has form 

in'^^K.[-r?vr"nr,). (37) 

While changing reference frame to K' , old wave function is written aJ^ 

1'(r,f) = (p(r',r')e''«"-'''\ (38) 

where (p denote wave function in new reference frame. Fractional derivative of product is 
revealed by generalized Leibniz rule: 

('P^'-^) = E w ^^^T^l^T 7-.DT"^Dyf. (39) 

„=or(a-n + i)r(« + i) ^ ^ ^ 

After substitution eqn.([38]l in ( [37| with taking into account eqn. ( [39| one can derive 
{ihd,i(p — Jicpdfif — ihvVyi^ + hvcpV/f) e'-^ = 

^.(-'■y^\|. r(a-' + orUi) °°'>°--»"- 



For Galilean invariance of eqn. ( [40] l one has to keep only time derivative and spatial 
fractional derivative of order a. Let us equate to zero coefficient near the rest deriva- 
tives of wave function (p{r,t). Then we can obtain that v — 0,(9,// = 0, V^' f — 0, that is 
f{r',t') = const. It leads to 



December 15, 2010 1:30 WSPC/INSTRUCTION FILE Paper' spmcs 



10 

where A is some constant. Thus, fractional Schrodinger equation is Galilean noninvariant, 
that is it changes its form with changing inertial reference frame. As a consequence one 
needs to choose the especial reference frame where, for example, nanoporous media is in a 
rest as well as to introduce additional potential into two-fluid hydrodynamic model. 

6. Conclusion 

It was shown that fractionalized set of hydrodynamical equations with taking into account 
nonextensivity of He-II inside nanopores leads to coupling between the first and the sec- 
ond sounds that appears even in the absence of viscous friction, which is undoubtedly the 
effect of fractional dimensionality of nanopore space and nonextensive nature of helium 
droplets. It was proposed that for better description of behavior of superfluid in nanoporous 
media with complex fractal structure one can use fractional Schrodinger equation. It leads 
to Galilean noninvariance and as a consequence one needs to choose the especial reference 
frame where, for example, nanoporous media is in a rest as well as to improve two-fluid 
hydrodynamic model by introducing there additional potential. 
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